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1  Introduction 


In  recent  years,  several  variations  of  the  Discontinuous  Galerkin  Methods 
(DGM),  for  second  order  elliptic  boundary  value  problems  have  been  pro¬ 
posed  which  exhibit  special  convergence,  conservation  and  local  approx¬ 
imation  properties  attractive  for  parallel  adaptive  hp-a p p ro x i m a  ti  on s .  An 
account  of  several  types  of  DGM's  can  be  found  in  the  book  edited  by  Cock- 
burn,  Karniadakis  and  Shu  [4], 

In  previous  work  [9],  we  introduced  a  new  stabilized  DGM  formula¬ 
tion.  Existence  and  uniqueness  of  stable  solutions  were  established  and  a 
priori  convergence  estimates  on  the  approximation  error  were  derived  for 
the  case  of  a  reaction-diffusion  type  model  problem. 

We  continue  our  error  analysis  of  this  stabilized  DGM  here  by  deriv¬ 
ing  global  implicit  a  posteriori  estimates  of  the  approximation  error.  Results 
of  our  previous  analysis  enable  us  to  prove  equivalence  between  the  error 
norm  and  the  norm  of  a  residual  functional  that  characterizes  the  accu¬ 
racy  of  the  approximate  solution.  Hence,  upon  estimating  the  norm  of  this 
residual,  we  indirectly  obtain  an  estimate  of  the  error.  To  verify  the  effi¬ 
ciency  of  our  estimates,  we  present  ID  and  2D  numerical  experiments  on  a 
reaction-diffusion  type  model  problem. 

In  Chapter  2  we  introduce  the  notations,  the  model  problem,  and  the 
stabilized  DGM  formulation.  Error  estimators  are  derived  in  Chapter  3 
and  numerical  verifications  are  presented  in  Chapter  4.  Finally,  concluding 
remarks  are  summarized  in  Chapter  5. 
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2  Model  Problem  and  Notations 


2.1  The  Reaction-Diffusion  Equation 

Let  C  R3  be  a  bounded  open  domain  with  Lipschitz  boundary  Oil  and 
let  {Ph}  be  a  family  of  regular  partitions  of  0  into  open  elements  K,  with 
diameters  hx,  such  that  (see  Fig.  2.1): 


=  int  |  |^J  K 

\Kevh 

The  maximum  diameter  in  the  partition  is  denoted  h.  The  set  of  all  edges 
of  the  partition  V},  is  given  by  £/,  =  { 7/s } .  k  =  1 .... .  iVedge/  where  iVedge  rep¬ 
resents  the  number  of  edges  in  the  partition  Vh •  The  interior  interface  Tint 
is  then  defined  as  the  union  of  all  common  edges  shared  by  elements  of  the 
partition  Vh'- 

iVedge 

hint  =  \Jik\dn. 

k= 1 

The  definition  of  the  unit  normal  vector  n  on  each  7^  is  related  to  the  num¬ 
bering  of  the  elements  in  the  partition,  such  that  n  is  defined  outward  with 
respect  to  the  element  with  the  highest  index  number  (see  Fig.  2.1).  The 
normal  vector  //  is  defined  outward  to  each  element  individually.  Within 
this  setting,  the  following  reaction-diffusion  problem  is  considered: 


—A u  +  u  =  /,  in  fi, 
u  =  0,  on  dU, 


(2.1) 


where  /  is  a  real-valued  function  in  L2(fl).  For  the  sake  of  clarity  in  the 
notation,  the  jump  and  average  operators  on  each  7^  E  Tint  are  respectively 
defined  as  (see  Fig  2.1): 


M  V\jkC3Ki  V'-tkC.»Kj 
(' V )  _  +  V  -fkC  t>Kj  ) 


7 k  =  int (dKi  n  dKj),  i  >  j. 


(2.2) 
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Figure  2.1:  Geometrical  definitions. 


2.2  The  Variational  Formulation 


In  this  section,  we  present  our  weak  discontinuous  formulation  of  the  model 
problem  (2.1).  First,  we  introduce  the  following  broken  space: 

M(Vh)  =  ([vtL2(n):  vlK  €H(A,K),VK€Vh,[Vvn]  eL2^)}, 

where 


H(A,K)  =  jt>  £L2(iL)  :  V- Vu  £T2(iL)J  C  H^K). 

Notice  that  v  £  H(A,K)  implies  S7v  ■  p  £  H~1!2{dK)  (see  [2,  7]).  The  norm 
I'll  on  M(Vh),  is  defined  as: 

Nil2  =  ^2  jll^llHhx)  +  GJ llVu  '  /illH-1/2(ax)| 

1  J  (2.3) 

+  H[V^-n]lll2(rint)- 

where  the  parameters  v,  A,  6,  and  (  are  all  positive  numbers.  At  this  stage 
of  our  development,  p  is  merely  a  non-zero  positive  integer  arbitrarily  as¬ 
signed  to  each  element.  Later,  this  parameter  will  be  taken  as  the  polyno¬ 
mial  degree  pk  of  the  approximations  on  each  element  K,  or  globally,  as 
the  minimum  value  of  all  px' s  in  the  partition  P/,  (see  Section  2.3).  Finally 
the  norms  in  (2.3)  are  defined  as: 


\\S\\H-1/2(dK) 

\\f\\H1/2(dK) 


sup 

ipeH1/2(dK) 


1(0!  <P)-l/2xl/2,dK\ 
\\iP\\H1/2{dK) 


inf  JHIffiutt, 

weHl{K)  v  J 

10W—1P 


(2.4) 
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where  (v)-i/2xi/2 ,8K  is  the  duality  pairing  in  H  1!2(dK)  x  H1!2( dK )  and 
where  70  denotes  the  trace  operator: 

70  :  H\K)  ~^H^2(dK). 


The  choice  for  the  space  of  test  functions  V  for  the  weak  formulation  is 
the  completion  of  M.(Vh )  with  respect  to  the  norm  |||•|||.  The  discontinuous 
variational  formulation  can  then  be  stated  as  follows: 


Find  w  €  V  : 

B(w,v)=C(v),  V'uGk 


(2.5) 


where  the  bilinear  form  B  :  V  x  V 
defined  as  follows: 


and  linear  form  C  :  V 


are 


B(w,v)  =  ^  <  /  {Vw  ■  Vv  +  wv}  dx 

Tftz'ry,  v  J  K 


Kevh 


{r(Vw  •  fi)  —  (Vv  ■  n)wj  ds 


Kevh 


} 


£(v)  =  f  fvdx. 

Jn 


/  {(v)  [Vw  •  n]  —  (w)  [Vv  •  n]  }  ds 

f  /jA 

+  L  w  ■  n]  [Vu  •  n]ds, 

J  Tint  P 


(2.6) 


The  following  notation  is  used  to  denote  the  duality  pairing  in  iT  1/2  (dK )  x 
Hl!2(dK)\ 

(Vw  ■  [i)v  ds  =  { Vw  ■  jl.  v)  H~1l2(8K)xH1l2(dK)- 

As  noted  in  [9],  this  weak  formulation  is  closely  related  to  the  DGM  for¬ 
mulation  by  Oden,  Babuska  and  Baumann  [6].  A  distinctive  difference  is 
the  addition  of  the  stabilization  term  on  the  jumps  of  the  normal  fluxes 
across  the  element  interfaces.  This  stabilization  allows  us  to  prove  the  well- 
posedness  of  the  weak  problem  while  having  the  advantage  of  not  pollut¬ 
ing  the  local  conservation  property,  which  is  one  of  the  appealing  properties 
of  DGM's.  Indeed,  by  taking  v  =  1  in  (2.6)  on  an  element  K  €  Vh,  we  obtain: 


/  wdx  —  /  (Vw  ■  p)ds  =  /  fdx, 

Jk  J  dK  Jk 
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and  local  conservation  is  satisfied  in  terms  of  the  average  fluxes  on  the 
element  edges.  The  well-posedness  of  the  problem  (2.5)  is  established  in  [9]. 


2.3  The  Discrete  Problem 

Let  {F^}  be  a  family  of  invertible  maps  defined  for  a  regular  partition  Vf, 
such  that  every  element  K  EVh  is  the  image  of  F  /x  acting  on  a  master  ele¬ 
ment  K,  i.e. 

Fir  :  K  — >  K,  x  =  F*(x).  (2.7) 

In  the  computational  model,  a  finite  dimensional  space  of  real-valued  piece- 
wise  polynomial  functions  of  degree  less  than  px  is  introduced  as: 

Vhp  =  £i2(H)  :  v]k  =DoF-\  v  E  PPK{K),  V  K  EVh}.  (2.8) 

We  note  that  Vhp  is  a  subspace  of  V.  Using  the  Galerkin  method,  an  approx¬ 
imation  Uh  of  u  is  sought  as  the  solution  of  the  following  discrete  problem: 
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3  A  Posteriori  Error  Estimation 


In  this  section,  we  derive  implicit  a  posteriori  estimators  of  the  approxima¬ 
tion  error.  These  are  global  in  the  sense  that  they  estimate  the  error  in  a 
global  norm  defined  on  the  whole  domain  0. 

We  present  two  types  of  estimators.  Section  3.3  introduces  the  estimator 
Xh  which  is  the  solution  of  a  global  problem.  This  type  of  estimator  will  be 
referred  to  as  a  coupled  estimator,  since  it  requires  solving  a  global  system 
of  linear  equations.  On  the  other  hand,  the  estimator  proposed  in  Sec¬ 
tion  3.4  is  obtained  by  solving  a  set  of  local  problems  and  is  referred  to  as  a 
decoupled- type  estimator. 


3.1  The  Error  Problem 

Let  and  u/ }  e  Vhp  be  the  exact  and  approximate  solutions  to  (2.5)  and 

(2.9),  respectively.  Then  by  using  the  linearity  of  the  bilinear  form  B( ■,  •),  it 
follows  easily  that  the  approximation  error  e  =  u  —  is  governed  by: 


B(e,v)  =  £(v)  —  B(v,h,v),  Vw6T 

V - v - ' 

nh{v) 


(3.1) 


where  T&  :  V  — >  R  is  the  residual  functional.  Note  that  due  to  (2.9)  the 
residual  satisfies  the  following  orthogonality  property  on  Vhp: 

Kh{v)  =  B(e,vh)  =  0,  V  vh  £  Vhp.  (3.2) 

Lemma  3.1 .1  The  residual  functional  is  continuous  on  V,  i.e.: 

300:  \Uh(v)\<C\\\v\l  VueV. 

Proof:  This  lemma  follows  trivially  from  the  continuity  properties  of  the 
variational  forms  £(•)  and  B(-,  •)  on  V  (see  [9]).  ■ 
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3.2  Equivalence  between  the  Residual  and  the  Error 


Since  lZh  is  continuous  on  V  with  respect  to  the  norm  Ill-Ill,  we  define  the 
norm  of  the  residual  (dual  norm)  as: 


sup 

oev\{0} 


\K\v)\ 

Nil 


Lemma  3.2.1  The  norm  of  the  error  and  the  norm  of  the  residual  functional  are 
equivalent,  i.e.  there  exist  two  positive  constants  C\  and  C-2  zvhich  depend  on  the 
mesh  parameters  h  and  p  such  that: 

C\  1  e  1  <  |||ftft|||,  <C2|||e|||. 


Proof:  Firstly,  by  recalling  the  Inf-Sup  property  of  &(-,■)  in  Theorem  3.3.2 
in  [9],  we  know  that  there  exists  7 (h,p)  >  0  such  that: 


7  (h,p) 1 e 1  < 


sup 

^ev\{o} 


\B{e,v)\ 


Using  (3.1)  and  the  definition  of  the  norm  |||7\F|||*,  we  obtain: 

7 (h,p)  1  e  1  <  sup  ^  ^  =  11^ I*. 
^ev\{o)  IF  III 


Conversely,  by  recalling  the  continuity  property  of  the  bilinear  form,  as 
stated  in  Theorem  D.l,  we  get: 

lll^lll*  =  sup  <  M(h,p) |||e|||, 

^ev\{o)  IF  III 


where  M(h,p)  =  max 


In  principle,  this  lemma  justifies  estimation  of  the  error  norm  indirectly  by 
estimating  |||7Nll*-  However,  it  is  essential  that  the  constants  C\  (h,p)  and 
C2  (h,  p)  do  not  degenerate  or  diverge  as  h  — >  0  and  p  — >  00.  This  implies  that 
both  the  Inf-Sup  and  continuity  constants  of  £>(■,•)  have  to  be  independent 
of  h  and  p.  According  to  Corollary  3.1  in  [9]  and  Theorem  D.l,  both  coeffi¬ 
cients  are  constant  if  A  =  u  =  (  =  9  =  0.  In  the  following,  these  parameters 
are  set  equal  to  zero  so  that  the  norm  |||- 1||  is  now  defined  as: 

Nil2  =  Y,  {llullHM^)  +  llVw'/illH^/2(ax)}+CTll[v«-n]llia(i«iIlt)-  (3-3) 

Kerh 
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3.3  A  Coupled  Error  Estimator 


A  natural  approach  for  error  estimation  could  be  to  approximate  (3.1)  using 
a  space  Whp  such  that  Vhp  C  Whp  C  V.  A  possible  candidate  for  Whp  is  the 
finite  element  space  of  functions  Whp  of  piecewise  continuous  polynomials 
that  are  of  degree  one  higher  than  the  functions  in  Vhp: 

Whp  =  {v  G  L2{n)  :  v\K  ='0oF-1,  vePPK+1{K)  V  A  G  Vh}  . 

In  this  case,  an  approximation  fh  G  Whp  of  the  error  is  obtained  as  the 
solution  to  the  error  problem  (3.1)  in  the  space  Whp,  i.e.: 


B(fh,vh)  =  llh(vh),  Vt)/,  G  Whp 


(3.4) 


and  a  global  estimate  of  the  error  is  given  by  the  norm  of  fh,  i-£-  |||  <f>h  III  ~  |||e 
as  indicated  in  the  following  lemma: 


Lemma  3.3.1  Let  G  Whp  be  the  unique  solution  to  (3.4),  G  Vhp  the  solu¬ 
tion  to  (2.9),  u  G  V  the  solution  to  (2.5),  and  let  u\K  G  H2(K)for  every  K  G  Vh- 
Then  if  the  family  of  mappings  (2.7)  is  affine  and  invertible,  there  exists  a  constant 
C(a )  independent  of  h  and  p,  such  that: 


e-fhl  <  C{a) 


hu-3/2 
pr~  7/2 


p  >  1,  rK  >  2, 


inhere  u  =  min(u  +  1,  r)  and  r  =  min  (r k)- 

Kerhy 

Proof:  See  Appendix  A.  ■ 

For  p  >  1,  I  ^  I  provides  an  asymptotic  estimate  of  the  error  norm.  The 
predicted  convergence  rates  are  suboptimal,  but  convergence  to  the  error 
norm  is  guaranteed  even  for  p  =  1.  Note  that  in  this  case  the  error  itself 
does  not  necessarily  converge  as  h  tends  to  zero  (see  Figure  5.1a  in  [9]). 

However,  although  this  approach  generally  provides  for  accurate  error 
estimates,  solving  the  resulting  system  of  equations  is  far  more  expensive 
than  solving  the  original  problem  for  u^.  Moreover,  for  linear  equations, 
this  approach  is  equivalent  to  computing  a  new  finite  element  solution 
Uh  £  Whp  of  (2.5)  and  taking  the  difference  Uh  ~  uh-  In  order  to  reduce 
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the  computational  cost,  it  is  desirable  in  a  posteriori  error  estimation  to  set 
up  local  problems  defined  on  each  element  or  patch  of  elements.  The  DG 
formulation  seems  well  suited  to  do  so.  Indeed,  let  the  test  function  v  be 
zero  everywhere  except  in  element  K.  The  equation  of  problem  (3.1)  thus 
reduces  to 


/  {Ve  •  Vv  +  ev}dx  —  /  {«(Ve  •  p)  —  (Vv  ■  p)e}ds 

JK  JdK 

+  [  {^[Ve-n]  -  (e)(Vu-n)}ds  (3.5) 

JdK/\dn  1 

+  /  cr[Ve  •  n]  (Vv  ■  n)ds  =  Hh x{vh) 

JdK/\an 

and  is  thus  defined  locally  except  for  the  coupling  terms  [Ve  •  n]  and  (e). 

As  an  intermediate  step  toward  full  decoupling  of  the  system,  we  intro¬ 
duce  the  function  x  €  V  that  is  governed  by: 


<A{x,v)  =  Tlh{v),  VveV 


(3.6) 


where  the  bilinear  form  A  :  V  x  V 


is  now  defined  as: 


A(w,v)=  <  /  {Vw  •  Vu  +  wv}  dx  > 

Kc-V,  K  J 

{v(Vw  ■  p)  —  (Vv  •  p)w}  ds 

+  /  dv  w  ■  n]  [Vv  •  n]  ds. 


Kevh 

E 


Kevh 


Comparison  of  A{-,  ■)  with  the  bilinear  form  B(-,  ■)  reveals  that  the  terms  in¬ 
volving  (v)  [Vw  •  n]  and  (■ w )  [Vv  •  n]  have  been  eliminated  from  the  formu¬ 
lation.  We  show  in  the  following  that  A(-,  ■)  satisfies  the  Inf-Sup  condition 
and  is  continuous. 


Lemma  3.3.2  (Inf-Sup  condition)  There  exists  7.4  >  0,  independent 
such  that: 


sup 

^ev/{o) 


\A{u,v)\ 


>  7.4  INI 


V«ey. 


of  h  and  p 


Proof:  See  Appendix  B. 
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Lemma  3.3.3  The  bilinear  form  A{-.  •)  is  continuous  on  V  x  V,  i.e.: 

3M_ 4  >  0  :  |*4(tt,v)|  <  M ^  |||rt|||  |||v|||,  V  u,v  E  V, 

where  M ^  =  2. 

Proof:  By  recalling  the  definition  of  .A(v)  and  applying  the  Schwarz  In¬ 
equality,  we  get: 

A{u,v)  <  E  iwiff  hiCjIklliJhK)  +  ||V«  •  h\\H-1/2(dK)\\v\\H1/2(dK) 
kevA 

+  11  Vv  •  h\\H-1/2(dK)Wu\\H1/2(dK)\  ^-7) 


+u  II  [Vu  ■  n]  I  ,,2(l  int)  ||  [Vv  •  n]  ||L2 (rint). 


By  definition  of  the  II 1  2  ( <9iL )  -nor  m,  we  know  that: 

\\V\\H1/2(dK)  —  K’l  // 1  (  K' } )  v  ^  Hl(K). 

Substitution  of  this  inequality  into  (3.7)  finishes  the  proof.  ■ 


Theorem  3.3.1  Let  x  €  V  be  the  unique  solution  to  (3.6).  Then  the  norm  |||x||| 
and  the  norm  of  the  residual  functional  |||7?/t|||*  are  equivalent: 


7.4 Ml  SIR'- III. 


Proof:  By  substituting  (3.6)  into  the  definition  of  |||7?/t|||*,  one  obtains: 

I*‘I.=  Sup 

vEV/{0}  IIMII 


(3.8) 


The  assertion  is  then  established  by  applying  Lemmas  3.3.2  and  3.3.3.  ■ 

Hence,  by  estimating  ||y||  we  indirectly  estimate  |||7?/t|||*  and,  therefore,  the 
error  norm  itself.  The  coupled  error  estimator  r\c  is  then  defined  by  the 
approximate  solution  Xh  £  Whp  of  x  such  that: 


Vc  = 


(3.9) 


and 


A(Xh,vh)  =  IZh(vh),  V  vh  G  Whp 


(3.10) 


We  prove  in  the  following  that  r]c  provides  for  an  acceptable  error  estimate 
of  11x10  and,  a  fortiori,  of  the  error  norm  itself: 
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Lemma  3.3.4  Let  x  €  D  be  the  solution  to  (3.6),  x\k  be  in  H2{K)  for  all  K  e 
Vh,  and  Xh  £  Whp  be  the  unique  solution  to  (3.10).  Then  if  the  family  of  mappings 
(2.7)  is  affine  and  invertible,  there  exists  a  constant  C(a )  >  0  independent  ofh  and 
p,  such  that: 


z,/z  3/2 

|||v  vi.ll  ^  C(er\ 

1  ^2  \\x\\ 2Hr{Ky  P>l,rK>2, 

Kevh 

lllx  X/illl  —  c'(  )  pr-7/2  y 

inhere  u  =  min(n  +  1,  r)  and  r  =  min  (rx)- 

Kevh' 

Proof:  See  Appendix  C.  ■ 

Although  the  estimator  r/c  is  global,  it  can  still  be  decomposed  into  a  sum 
of  element-wise  contributions  r)c,x  such  that: 

vl,K  =  WXhW'ni (k)  +  W^Xh  '  T I \2h-  1/2 +  2a\\  i^Xh  ■  n]  Ilx,2(ajcnrint) 

so  that 

2  \  2 
he  =  2^  hc,K- 

Kevh 

The  contributions  rjc.x  can  serve  as  refinement  indicators  for  each  element 
in  an  mesh  adaptation  strategy. 


3.4  A  Decoupled  Error  Estimator 

In  this  section  we  seek  to  derive  an  estimator  which  involves  solving  local 
problems  only  on  each  element  of  the  mesh  in  order  to  reduce  the  compu¬ 
tational  cost.  Decoupling  the  global  equation  supposes  to  eliminate  the  in¬ 
teractions  between  neighboring  elements  through  the  element  boundaries. 
In  doing  so,  a  loss  of  accuracy  with  respect  to  the  estimator  r]c  is  therefore 
expected. 

To  establish  the  decoupling,  we  decompose  the  finite  dimensional  space 
Whp  into  the  family  of  local  spaces  {IT^}: 

W*g=  jv  =  t)o  F~\  V  €  Ppk+1(K)Y 

with  associated  norms: 

Mwk  =  Mm(K)  +  II  Vv  •  pW'^2 +  2a  ||  Vu  •  p\\2L2^dKr[Tiniy  (3.11) 
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Using  this  function  space  setting,  we  then  consider  the  following  local  prob¬ 
lem  for  every  K  £  V^: 


Find  V4  G  W«p  : 

VK(^hK,vh)  =  nhK(vh),  V  vh  G  1U^ 


(3.12) 


where  7ZhK{-)  denotes  the  restriction  of  7Zh(-)  to  K  and  V K ( • .  •)  is  the  bilin¬ 
ear  form  defined  on  Wfff  x  Wfff  such  that: 

T^K^K^h)  =J  \y^K^vh  +  ^Kvh\dx  +  (V^4  '  p)  •  fi)  ds 

+  2cr  /  (V^4  •  v){Vvh  ■  n)ds. 

JdKnTint 


Remark  3.4.1  There  exists  a  unique  solution  e  to  (3.12).  Indeed,  since 
is  finite  dimensional,  existence  and  uniqueness  of  ip \  follows  by  application 
of  the  Fredholm  alternative  and  by  noting  thatV( •,  •)  is  positive  definite  on  Wfff  x 
Wlfp. 

Once  again  we  can  consider  the  global  error  estimator  rpi  in  terms  of  local 
contributions  r),u<: 


so  that: 


Vd,K  =  HkWwk  =  ^VK(iphK,iphK) 

2  \  ^  2 

Vd=  2^  V d,K ■ 

Ken 


(3.13) 

(3.14) 


Theorem  3.4.1  Let  {ipx}  be  the  set  of  unique  solutions  to  (3.12)  and  let  %  be 
the  error  estimator  as  defined  in  (3.14).  Then  the  following  inequality  holds: 


%  <  IIIR'II. 


Proof:  Starting  from  the  definition  of  r/a  and  the  formulation  of  the  local 
problems  (3.12),  we  get 

rii=  £  Vk^k^k)  =  £  nhK{i>hK )  =  Kh{  £  v4)  O-is) 

Ken  Kevh  Ken 
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where  we  assume  that  each  local  function  is  extended  to  the  whole 
domain  such  that  riphK  =  0  outside  of  K.  For  simplicity  we  introduce  the 

global  function  iph  E  Whp  such  that  iph  =  'iphK .  It  follows  that: 

Kevh 

vl  =  Kh(i,h)  <  IK‘I.11^1. 

Let  Vh  be  an  arbitrary  function  of  Whp.  Now,  applying  the  Cauchy-Schwarz 
inequality,  we  observe  that: 

INI2  =  Y  {lM H1{K)  +  II' Vvh  ■  nfg-i /2{dK)  }  +  HI  [V^  •  n]  llL2(rtat) 
Kevh 

—  Yj  {ll^lltfhK)  +  W^Vh  ■  V\\2H-i/2(dK)} 

K£Th 

+  2crX/  {HV^^+llI2(^nrtat)  +  llV^^  llL2(aicnrint)} 

K£Vh 

—  X/  (KIlW)  +  IIV^  ■Mlll2^)  +  2^l|V^^|||2(aKnrint)} 

Kerh 

=  Yj  WVh\\wK- 

Kevh 

Replacing  Vf,  by  iph  in  the  above  inequality,  we  get: 

v2d  <  Ills'll* lll^lll  <  \\K%  I  Y  ll^ll wK  <  Vd\\\Kh\l 

V  K^h 

which  allows  to  conclude  the  proof.  ■ 
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4  Numerical  Results 


4.1  One-Dimensional  Experiments 


We  consider  the  one-dimensional  version  of  (2.1): 


d2u 

j  +  tt  =  l,  forO<x<l, 
u(0)  =u(l)  =  0. 

The  exact  solution  of  this  problem  is  given  by: 


(4.1) 


u(x)  =  1  — 


ex  +  e1-* 
1  +  e 


(4.2) 


In  the  following  experiments,  a  is  set  equal  to  one  and  the  quality  of  a  given 
error  estimator  r]  is  measured  in  terms  of  the  effectivity  index,  namely  the 
ratio  of  r/  and  |||e|||.  We  show  in  Fig.  4.1  the  effectivity  indices  for  rjc  and 
rjd  in  the  case  of  uniform  approximation  degrees  p  =  1,  2,  3  and  4,  and  for 
uniform  mesh  refinement.  For  comparison,  we  also  show  the  effectivity 
index  for  |||</>/J|,  that  is  the  ratio  |||^/i|||/|||e|||.  We  observe  that  the  estimate  r)c 
is  very  accurate  for  p>  2,  with  effectivity  indices  ranging  between  89  and 
100  percent.  For  p  =  1  however,  the  estimator  is  less  accurate  as  |||e|||  does 
not  necessarily  converge  in  this  case.  Nevertheless,  the  effectivity  index 
remains  bounded,  as  expected  since  |||e|||  and  |||y|||  are  equivalent  and  the 
fact  that  Xh  converges  to  y. 

The  results  for  the  decoupled  estimator  r/d  show  that  the  effectivity  in¬ 
dices  are  still  accurate  despite  the  decoupling.  They  are  within  1  or  2  per¬ 
cent  range  of  the  results  for  the  coupled  estimator  r]c  when  p  >  2.  However, 
for  p  =  1,  the  decoupled  estimator  is  not  as  accurate,  certainly  because  once 
again  the  norm  of  the  error  does  not  converge  to  zero  as  h  is  decreased. 


4.2  Two-Dimensional  Experiments 

For  our  2D  example  problem  we  consider  the  following  BVP,  given  on  the 
unit  square  =  (0,1)  x  (0,1)  with  prescribed  Dirichlet  boundary  condi- 
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Effectivity  Index 


2.5 


p  =  1  — i — 

p  =  2  — - x — 

p  =  3  * 


1.5 

1 


i — xx  r-xx-^t-x-x— x--x — x — I — -x 


0.5 


0  . . . . 

10  100  1000 

dofs 


Figure  4.1:  Effectivity  indices  versus  the  number  of  degrees  of  freedom  for  |||<^h||| 
(top),  for  the  coupled  estimator  rjc  (bottom  left)  and  the  decoupled 
estimator  %  (bottom  right)  -  ID  residts. 
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tions  on  dil: 


in  If, 

ye  [0,1], 


—A  u  +  u  =  0, 
u(0,y)  =  0,  u(l,  y)  =  0, 

ti(a;,0)=0,  u(x.  1)  =  -  sin(7nr)  sinh  \/l  +  7 r2, 
The  exact  solution  to  this  problem  is  given  by: 

u(x:y)  =  ^  sin(7nr)  sinh( \J\  +  7r2y ) . 


x  E  [0, 1]. 


(4.3) 


(4.4) 


As  before  the  parameters  A,  v,  (  and  9  were  set  equal  to  zero,  a  was  set 
equal  to  one,  and  computations  for  uniform  orders  of  approximation  p—  1, 
2,  3  and  4  were  performed  while  uniformly  refining  the  mesh.  However, 
the  calculation  of  the  H~1!2  (dK)- norm  in  two  dimensions  is  a  rather  com¬ 
plicated  task.  Therefore,  we  choose  here  to  approximate  the  norm  |||  •  |||  as 
follows: 

Nil2  ~  {ll^llffhic)  +  II }  +<JH['vV  n]  llz,2(rjnt)- 

KEVh 


Figure  4.2  shows  that  the  coupled  estimator  r\c  is  again  very  accurate  in 
two-dimensional  computations.  Its  effectivity  index  ranges  between  85 
and  110  percent  for  all  orders  of  approximation.  The  results  are  actually 
very  similar  to  the  results  obtained  for  |||^|||.  On  the  other  hand,  for  p  =  1, 
we  again  observe  an  overshoot  in  the  effectivity  index,  but  the  overshoot 
remains  bounded  and  is  far  less  important  than  is  the  case  for  the  one¬ 
dimensional  results.  Finally,  Figure  4.2  shows  a  dramatic  loss  in  accuracy 
for  the  uncoupled  estimator  r/(;  in  2D.  The  effectivity  indices  of  rj,i  ranges 
between  15  and  20  percent  for  p  =  1,  2  and  4  and  are  slightly  better  for 
p  =  3  with  values  converging  around  40  percent. 
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Figure  4.2:  Effectivity  indices  versus  the  number  of  degrees  of  freedom  for  |||<^h||| 
(top),  for  the  coupled  estimator  rjc  (bottom  left)  and  the  decoupled 
estimator  %  (bottom  right)  -  2D  residts. 
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5  Concluding  Remarks 


A  posteriori  global  error  estimates  are  derived  for  a  Stabilized  Discontin¬ 
uous  Galerkin  Method  (SDGM).  By  first  proving  an  equivalence  relation 
between  the  norm  of  the  residual  and  the  norm  of  the  error,  we  compute 
the  error  estimates  by  estimating  the  norm  of  the  residual  functional.  Two 
types  of  estimators  are  proposed:  1)  Coupled  estimator.  This  estimator  in¬ 
volves  the  solution  of  a  global  problem  on  the  entire  domain  due  to  the 
presence  of  coupling  terms  in  the  formulation.  It  is  obviously  expensive 
to  compute;  however  it  was  investigated  as  an  intermediate  step  between 
the  approximation  of  the  error  problem  and  a  fully  decoupled  error  estima¬ 
tor.  One-  and  two-dimensional  numerical  experiments  show  that  this  cou¬ 
pled  estimator  is  still  accurate  (although  certain  terms  have  been  eliminated 
from  the  error  equation)  and  generally  yields  effectivity  indices  ranging  be¬ 
tween  85  and  110  percent.  2)  Decoupled  estimator.  This  estimator  involves 
the  solution  of  a  set  of  local  problems.  One-dimensional  results  show  good 
accuracy,  with  values  of  the  effectivity  index  varying  between  80  and  100 
percent.  However,  in  two  dimensions  the  decoupling  leads  to  a  dramatic 
loss  in  accuracy  with  effectivity  indices  ranging  between  15  and  40  percent. 

More  numerical  experiments  would  be  necessary  to  understand  what  is 
happening  during  the  decoupling  procedure  in  two  dimensions.  However, 
we  believe  that  an  accurate  decoupled  error  estimator  can  be  derived  for 
the  present  version  of  the  Discontinuous  Galerkin  method.  An  alternative 
approach  would  be  to  set  up  the  local  problems  on  patches  of  elements,  in 
which  case  partial  exchange  of  information  would  be  allowed,  following  a 
subdomain-type  residual  method  proposed  for  continuous  finite  element 
approximations  (see  e.g.  [3,  5,  8]). 
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Appendix  A 

Proof  of  Convergence  Theorem  for  |||<4||| 


The  existence  and  uniqueness  of  4>h  E  Whp  to  (3.4)  follows  by  application 
of  the  Fredholm  alternative  and  remarking  that  the  bilinear  form  B(-, )  is 
positive  definite  on  Whp  x  Whp.  Now,  if  we  define  the  approximation  of  u 
in  the  augmented  finite  dimensional  space  Whp  as  u*h,  i.e.: 

B(u*h,vh)  =  £(vh),  VvhEWhp , 

and  substitute  this  equation  into  (3.4),  then  we  observe  by  the  linearity  of 
the  bilinear  form  that  fyh  =  u*h  —  u^.  Therefore, 

llle -41  =  14-411- 

We  now  recall  the  global  interpolation  operator  II/l},(-)  that  we  introduced 
in  our  previous  work  [9]: 

n /*  :  F  — >  Whp,  Uhp(v)  =  Y,  VT ■&(«,*),  (A.l) 

Kevh 

where  the  local  interpolants  4^(-)  are  defined  as  follows  [1]: 

7T^:  HrK{K)  — >  PPK+1(K), 


<,(«/.)=«/.,  Vvh€  PPK+1(K). 

Given  the  interpolant  II^p  (u),  we  split  a  —  u*h  such  that  u  —  u*h  =  r]  —  £,  where 
r]  =  u  —  II/j pu  and  £  =  u*h  —  H),pu.  Notice  that  £  E  Whp  and  that  the  interpo¬ 
lation  error  r]  E  H2(Vh)-  Thus,  by  using  the  triangle  inequality,  one  obtains: 

Il|e-4III<IIM|I  +  I4lll-  (A.2) 


By  applying  the  discrete  Inf-Sup  property  of  the  bilinear  form  in  The¬ 
orem  3.3.3  in  [9],  and  by  recalling  that  u  =  9  =  0,  we  know  there  exists  C  >  0 
such  that: 


141  SUP 

n  vhewhp/{ o) 


m,vh)\ 

lll^lll 


(A.3) 
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From  (3.1)  and  (3.4),  one  easily  observes  that  e  —  < $>h  satisfies  an  orthogonal¬ 
ity  property  on  Whp: 

B{e-4>h,vh)  =  0,  V  vhe  Whp. 


By  employing  this  property,  we  can  rewrite  (A.3)  as: 


III  £  III  <c  \  SUP 

n  vhewhp/{ o) 


\B(n,Vh)\ 

IK  III  ' 


If  we  apply  the  continuity  property  of  &(-,  •)  and  note  that  the  continuity 
coefficient  M  is  a  constant  due  to  our  choice  of  values  for  v  and  9  (see 
Section  3.2),  then  we  can  bound  |||£|||  by  the  interpolation  error  as  follows: 

Kill  <c(<k  III* 

Thus,  returning  to  (A.2)  we  conclude  that: 

lie  -  «  <  C(c)  £-  |h|||. 

By  recalling  the  interpolation  Theorem  4.2.2  in  [9]  and  remarking  that  v  = 
A  =  9  =  (  =  0  and  that  the  order  of  approximation  in  Whp  is  p  +  1,  we  get: 

h^-1/2  I - 

III7? Ill  —  r- 3/2  a  X/  WU\\2HrK{Ky  rK  —  2‘ 

7  V  Ken 

By  combining  the  last  two  results  we  establish  the  assertion  in  Lemma  3.3.1. 
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Appendix  B 

Inf-Sup  Property  of  Bilinear  Form  A(-,  •) 


The  proof  of  Lemma  3.3.2  is  analogous  to  the  proof  of  Theorem  3.3.2  in  [9]. 
Given  a  u  E  V,  then  we  can  construct  a  u  E,  such  that: 

u  =  u  +  f3  ^  zK,  (B.l) 

Kerh 


where  (3  E  R  and  the  local  functions  zk  £  17 1  [K)  PI  V  have  the  following 
properties: 

V zk  ■  H  =  Vti  •  /u,  on  <9iC, 

(z,v)  1,K=  (Vri  •  /j)  v  ds,  (B.2) 

JdK 

\\zK\\m{K)  —  l|V«  ■  ^\\H~1/2(dK)j 

where  (-,-)i :k  denotes  the  II1  (K)  inner  product.  For  more  details  on  the 
proofs  of  these  properties  and  the  construction  of  u,  the  reader  is  referred 
to  [2,  7,  9].  Now,  by  the  definition  of  the  supremum  we  get: 


sup 

vev/{o} 


\A(u,v)\ 


> 


\A(u,u)\ 


u 


(B.3) 


If  we  choose  (3  =  — 1/4,  for  instance,  then  by  simply  applying  the  definition 
of  A(-.  •),  we  can  expand  the  term  in  the  nominator  as  follows: 

A{u,u)  =  -  \{u,zk)i,k 

KEVh  1 

+  7/  {zk{Vu  •  /Li)  -  (VzK  ■  n)u}ds  \ 

4  JdK  J 

+  11  [Vu  •  n]  \\2L2iTint)  -  T  f  [V«  •  n]  [VzK  ■  n]ds. 


23 


Substitution  of  the  first  two  properties  in  (B.2),  gives: 


A(u,u)  =  ^  { 

KEVh  L 

+  |ll[V«-n]||i2(rtat) 


\m(K)  2 

|2 


-  ^(u,zk)i,k  +  ^ \\zk\ 


m(K) 


(B.4) 


By  recalling  Young's  Inequality: 


1 


(u,zr)i,k  <  2^u^h1(k)  +  2£ll Zr\\h1(k) 
and  choosing  e  =  2,  we  can  rewrite  (B.4)  as: 


,  e  >  0, 


A{u,u)  >  {^|||i(K)  +  ^||^|||i(X)]  +  ^||[Vri-n] 

IT  £1^0,  v  ' 


I  L2  (lint) 


KEVh 


By  substituting  the  third  property  in  (B.2)  into  this  expression,  one  obtains: 


A(u,u)  >  ^  |H||2. 

O 


(B.5) 


Conversely,  Lemma  3.3.1  in  [9],  yields  the  following  inequality  for  the  de¬ 
nominator  of  (B.3): 

1  ti  1  <  V2|«|.  (B.6) 

Finally,  substitution  of  (B.5)  and  (B.6)  into  (B.3)  yields: 

\A(u,v)\  ^  1  hi  in 

sup  — jr— TTj —  >  — -=  m  III. 

^ev/fo}  IF  III  8y2 
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Appendix  C 

Proof  of  Convergence  Theorem  for  |||x/J 


The  existence  and  uniqueness  of  Xh  £  Whp  to  (3.6)  follows  by  application 
of  the  Fredholm  alternative  and  remarking  that  the  bilinear  form  A( •,  •)  is 
positive  definite  on  Whp  x  Whp.  Now,  to  prove  the  lemma,  we  will  need  an 
discrete  Inf-Sup  condition  of  the  bilinear  form  A(-.  •)  on  Whp  x  Whp. 

Lemma  C.1  Let  {F^  }  define  a  family  of  affine  invertible  mappings.  Then  there 
exists  a  yh  =  7 h{h,p)  >  0,  such  that: 


\A{uhlvh)\ 

Kill 


>  7ft  IK  III, 


Vuh€Vhp/{  0}, 


and  7 h  =  C—t  >  0. 
jr 

Proof:  By  definition  of  the  supremum,  we  get: 

\A(uh,vh)\  ^  A(uh,uh) 

„6™?/f0)  III''* III  -  HI  ■ 

By  applying  Corollary  3.3.3  in  [9],  it  is  clear  that  there  exists  C  >  0,  such 
that 

A(uf „Uh)  >  cl  ^  (\\uh\\'2Hi(K)  +  72  W^Uh  '  nllff-i/2(ax) 

^K£Vh  V  1 

+a\\  [Vii/,  •  n]  ||i2|rjnt)  | * 

which  concludes  the  proof.  ■ 

Given  the  interpolant  Lippi')  in  (A.l),  we  split  x  —  Xh  such  that  X~Xh  =  V~ 
C  where  rj  =  x  ~  FI hpX  and  £  =  Xh  ~  Notice  that  £  e  Whp  and  that 
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the  interpolation  error  rj  E  II2  {Vh).  Thus,  by  using  the  triangle  inequality, 
one  obtains: 

Ix-xaI<M  +  KI-  (CD 

By  applying  the  discrete  Inf-Sup  property  in  Lemma  C.l,  we  know  there 
exists  C  >  0  such  that: 


III£III<c'y  sup 

n  vhewhp/{ 0} 


I -4(^)1 
IK  III  ' 


(C.2) 


From  (3.6)  and  (3.10),  one  observes  that  X  ~  Xh  satisfies  an  orthogonality 
property  on  Whp: 

A(x~Xh,vh)  =  0,  V  vh  E  Whp. 


By  employing  this  property,  we  can  rewrite  (C.2)  as: 


III£III<c'y  sup 

n  vhewhp/{ o} 


I  Av,vh)\ 
lll^lll 


If  we  apply  Lemma  3.3.3,  then  we  can  bound  |||£|||  by  the  interpolation  error 
as  follows: 

Kill  <  C(a)  £  Mil. 

Thus,  returning  to  (C.l)  we  conclude  that: 

IX-xdl<C(a)L|H|. 

By  recalling  the  interpolation  Theorem  4.2.2  in  [9]  and  remarking  that  v  = 
A  =  8  =  (  =  0  and  that  the  order  of  approximation  in  Whp  is  p  +  1,  we  get: 

up- !/2  I - 

Ml  <  C(a)—^  Hxll %rK(K),  rK  >  2. 

p  V  Kevh 


By  combining  the  last  two  results  we  establish  the  assertion  of  Lemma  3.3.4. 
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Appendix  D 

Continuity  Property  of  Bilinear  Form  B(-,  •) 


Theorem  D.1  Let  £>(., .)  be  the  bilinear  form  as  defined  in  (2.6).  Ifc r  >  0,  then: 


3M  >  0:  \B(u,v)\  <  M|||tt|||  |||v|||  Vu,vEV 


where  M  =  max  <3, 


tpO  I  PC  j 


hv''  V  2ahx  y 


Proof:  From  the  definition  of  £?(., .)  and  observing  that  for  u,  v  EV  we  can 
write: 

f  (v)  [Vtt  •  n]ds  =  ^  f  v  [Vtt  •  n]  ds,  (D.l) 

*  dKDTint 


KEVh' 


which  yields: 

B(u,v)  =  {Vrt  •  Vv  +  uv}dx 

■  /  v{S7u-n)ds+  /  u(Vv-ii)ds 
JdK  JdK 


KEVh 


1 

+  2 


j  ti[Vw  •  n]  ds  —  ^ 

/  -u[Vu  •  n]  ds  > 

JdKnTint  2  • 

JdKnTint  J 

hx  f 

+  <7—  /  [Vtt  •  n]  [Vv  •  n]  ds. 

P  J  Tint 


By  applying  the  Schwarz  Inequality  and  by  using  the  definition  of  the 
H 1/2  (dK )  norm  (2.4),  we  can  bound  the  above  as  follows: 
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